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ABSTRACT 


A method is presented for determining the relationships 
between the costs and technical performance of environmental 
vapor power systems in a manner which permits fundamental 
design specifications to be made optimally with respect To 
overall system lifetime costs. Means of applying optimization 
techniques for large scale systems to the thermoeconomic 
analysis of environmental vapor power systems are described 
and demonstrated with a simplified sample model. A sequential 
unconstrained minimization algorithm is employed for overall 


system design optimization. 
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A. BACKGROUND 

The motivaticn for developing new energy technology, long 
forecast by such researchers as Putnam [1], has now become so 
widely understood as to recuire no elaboration here. Research 
and development efforts are proceeding on a broad front in 
search of alternatives to the conventional non-renewable fossil 
fuels and potentially hazardous fission processes. 

One class of proposals seeks to extract useful energy 
from sources existing in nature. With the exception of geo- 
thermal energy, virtually all of these rely ultimately on some 
phenomenon associated with receipt by the earth of energy 
radiated from the sun. A sub-class of these "environmental" 
energy systems utilize vapor cycles in mechanizing the 
conversion from the diffuse heat sources found in nature into 
the more concentrated and transportable energy forms required 


in many applicaticns. 
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Figure 1 diagrams the fundamental concept upon which these 
environmental vapor power systems depend. Heat is extracted 
from some natural source of elevated temperature (geothermal 
wells, direct solar collectors, hot seawater, etc.) and 
transferred to a working fluid. Devices suitable to the 
application convert the thermodynamically available portion of 
the available heat into other forms of energy, and the 
unavailable energy 1s rejected to a thermal sink. Both open 
and closed vapor cycles are possible [2], [3] and the products 
of conversion can take many forms, including electricity and 
fuels such as hydrogen, methanol, and ammonia [4]. 

One widely used vapor cycle is the closed Rankine cycle 
[5], shown schematically in Figure 2. 

In addition to the conservation of chemical fuels, 
certain of the proposed methods of harnessing energy hold 
promise of significant additional advantages. It is expected 
that their effects on the environment will be relatively benign 
[6], particularly in terms of atmospheric and thermal pollution. 
They would cause no addition to the total heat burden at the 
earth's atmosphere and rely on sources which are continuously 
renewed by natural processes [7]. 

There 1s another class of vapor power cycles which, 
although not always exploiting environmental energy sources, 
shares enough of the operating characteristics of those that 
do to warrant mention as a group amenable to the type of analy- 
Sis discussed in this paper. These are the "bottoming" cycles 


for extracting power from the still energetic discharges of 
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geothermal, nuclear, and chemically fueled plants. Figure 3 
shows two general types of these "waste heat"cycles. 

The essential feature which differentiates both the 
environmental and waste heat vapor systems from conventional 
ones is the relatively low thermal potential within which 
they operate. A consequence of this characteristic is that 
the size scale of all the cycle components is increased in 
comparison with conventional plants. Heat exchange surface 
areas must be enlarged for sufficient heat to be transferred 
through small driving potentials, and with less energy avail- 
able from each unit of fluid circulated, a far greater volume 
rate of working fluid must be cycled. Pumps, pipes, and 
conversion devices such as turbines all grow in size as the 
temperature difference between the source and sink is reduced 
while the energy product is held constant. 

Viewed fundamentally, environmental power systems employ 
technology which has been available for many years. Many 
concepts have been tested with working models or demonstration 
plants, and some are employed presently on a smali scale. 
Although significant technical problems arise in connection 
with specific applications, these do not appear to be permanent 
obstacles. Net energy asessments appear favorable and 
questions of material availability and local adverse environ- 
mental effects seem amenable to solution [6]. 

The primary question which will determine when environ- 
mental energy sources will be exploited on a scale large 


enough to significantly affect the energy market is that of 
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system economics. Although some researchers predict plant 
costs which are currently competitive with conventional 
methods, [9] uncertainties arising from the lack of operating 
experience weaken the claims of these proponents. Until 
economic viability can be conclusively shown, risk aversity 
will act as a strong deterrent to attracting the very large 
amounts of venture capital required. With the private sector 
presently unconvinced, the federal government is undertaking 
the expenditures required for research and development efforts 


eo]. 


B. | THERMOECONOMICS 

Typically, system design and cost studies are conducted 
in a two-step or, at best, iterative process. Designers 
assemble specifications based on technically achievable and 
desirable functional characteristics. They are, of course, 
guided in their design decisions by some measure of intuition 
as to the economic impacts, usually based on prior experience 
with similar programs. The degree of detail in the initial 
specifications presented, in fact, often reflects the confi- 
dence held by the engineers in their economic appraisals. The 
system and component specifications are then subjected to cost 
analysis, prime cost factors are identified, and technical- 
economic trade-offs are suggested. 

The probleins arising from this partial separation of the 
design and costing steps are more or less severe according to 
the application. There is a fundamental difficulty in 


communicating the two groups' understandings in a meaningful 


ie 





way, a difficulty which increases as the novelty of the design 
Situation and hence number of unconstrained design choices 
increases. In their quest for a sophisticated design product, 
engineers may design around some apparently desirable 
Parameter, such as a high heat transfer coefficient. Cost 
analysts may take this figure as fixed and address themselves 
to questions of material selection, maintainability, or 
manufacturing tolerances without recognizing that adjusting 
the heat transfer coefficient itself could produce the most 
rewarding cost effects. 

This type of difficulty is most severe when little ex- 
perience is available to guide the engineer's fundamental 
choices, as is true in the case of environmental vapor power 
systems. The small available thermal potential drastically 
distinguishes these systems from their high temperature 
counterparts. Ina fuel-fired generating plant, for example, 
the power required to pump the working fluid can be neglected 
in a first approximation, and a variation of 1°F temperature 
difference across a boiler tube is insignificant. As is 
demonstrated later in this study, such considerations can have 
a profound influence on the overall economics of an environ- 
memtal plant. 

To some extent, the engineer's problem can be viewed as 
being where to start. Recognizing that pumping power 1s going 
to be substantial, he might decide to assign 10 or 20 percent 
of the plant's overall output to pumping requirements and 


build much of the rest of the design about this choice. Or 


1s 





he might choose to utilize 30 percent of the total temperature 
difference for heat transfer, leaving the remaining 70 percent 
available for enthalpy drop across the turbine. He might 
establish a dimensional constraint, based on nothing much more 
than the feeling that a 100 foot diameter pipe is a very big 
pipe. 

Unfortunatly, all these basic choices involve performance 
and cost tradeoffs. If flow rates are increased to enhance 
heat transfer, drag coefficients increase as well. How much 
improvement in heat transfer is worth how large an increase in 
pump head, and hence pump work? Pump work is also influenced 
by heat exchanger tube diameter, Spacing, and surface charac- 
teristics, which also affect space and nee requirements. 
How much should one be willing to pay to reduce fouling heat 
resistance? If heat exchange 1s dominated by fouling 
resistance, is it worth the extra temperature drop necessary 
to shift to a different boiling regime? Unless the cost 
analyst is knowledgeable about the thermodynamic consequences 
of costing factors he is in as poor a position as the engineer 
to make the tradeoffs in dollars per millimeter of fouling 


organisms. 


OF OE weeny E 

What is needed is an analytical method whereby overall 
economic effects may be integrated into engineering design in 
such a way that the designer's intuition may be enhanced in 
trading off the costs and benefits of parameter selection at 


the margin. A means is required for mapping the large number 





of interrelated engineering variables into their individual 
and collective effects in the marketplace, where the ultimate 
design appraisal will take place. 

The research reported on in this paper is intended to 
develop and evaluate one method of integrating marginal 
cost/benefit analysis into engineering design and to show 
the kinds of information which could thus be gained. In this 
initial investigation, no effort has been made to apply the 
method to any particular practical design problem or to 
produce analytical insights into existing systems. The 
intent has been to show how thermoeconomic analysis can be 
performed and what value it can have when applied to a specific 


real case. 
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If. THERMOECONOMIC ANALYSIS 


A. CONCEPT 
Profit is the difference between benefits and costs, 
both broadly considered. When these can be related over a 


~ 


common set of decision variables, X, one may write 
w(X) = BCX) - CCX) 


with B(X) representing the sum of all benefits, and C(X) the 


Sumeor all costs: 


aS N s 
B(X) = & B.(X) 
a ab 
a N Pe 
CO) mens ae. (X) 
a 


If all the relevant B, and C; can be defined functionally 


~e 


over X, performing 


maximize: 7 


subject to: a required level of performance (A) 


would produce the desired optimization. 


B. Gilt REDUCTION 

Attempting a global optimization directly with all possible 
costs and benefits considered, although theoretically possible, 
encounters many practical difficulties [11]. It is possible, 
however, to achieve considerable reduction of the problem 


without sacrificing many of the benefits of the analysis. 
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First, although the impacts of general (and difficult to 
quantify) externalities, such as independence from foreign 
control of energy sources, are important and should not be 
excluded from the final analysis, many interior decisions 
suffer not at all from excluding externalities such as these 
from most of the study. Many other factors are not related 
to the decision variables (Xs) and therefore do not affect 
marginal design choices. For example, personnel training 
expenses are not close functions of tube diameter. Since 


miemsolution to 


Maxtimezesm = BCX) = CCX) = D 


subject to g(X) = 0 


ad 


where D is constant with respect to X is identical to the 
Bomuti1on of 
maximize mt = BCX) = CCX) 


subject to - g(X) = 0 


any factor which acts only as an additive constant may be 
excluded from the analysis without affecting the results. 
Even with invariants over the decision variables ex- 
cluded, there are other serious impediments to seeking global 
solutions to (A). Convexity of the optimization problem is 


not assured by the physical relations modeled.* When ali 


Thy convex optimization problem is defined as one with a 
convex objective function, to be minimized,concave ? inequal- 
ity constraints, and linear equality constraints [12]. The 
conditions on the constraints assure that the feasible region 
Mmemcmeconves Gel, 1.c¢., for every X%, 0 < A < 1, and any two 
points X,, X, €T, a convex set, [AX + (1-A)X, 1 er, 





decision variables are considered at once in a global assault 
it is increasingly difficult to test for uniqueness of the 
solution. Secondly, since design variables in one system 
component often are only distantly related to those in another, 
insights are obscured when they are varied simultaneously 
within one code. Thirdly, the model can never be exact. It 
is important for the designer to keep track of the effects 
of his.modeling choices in detail. This is more easily achieved 
by putting the pieces together sequentially than all at once. 
Finally, the designer often has adequate information available 
to intelligently fix some of the variables. It is unnecessary 
to complicate the analysis by including as free variables 
factors which are closely constrained by other considerations. 
For these reasons, it appeared desirable to fellow the 
usual procedure for the optimization of large scale systems 
by decomposing the problem into coherent interrelated zones 
and achieving global optimality through one of the available 
zone coordination methods. The next section contains an 


outline of the general procedure. 


c. PROBLEM COORDINATION 

The general theory for optimizing large scale systems 
through coordination of smaller subsystems can be found in 
references such as Wismer [11]. The following discussion of 
the two basic approaches is greatly particularized in that 
the terminology and composition of the examples reflect the 
structure of the sample analysis which is presented in section 


Ho]. 
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The first approach, called the model coordination method, 
can be understood through consideration of the decomposed 
system shown in Figure 4. 

Define . = (Ti,2TQ) as a vector of coordinating variables 
and X, and - as vectors of design variables in zones 1 and 2. 
Then construct the zone subproblems: 

minimize: £.(X,5y) 

X- sy sus enellsc2 


subject to: g.(X.,y) > 0 


tl 
© 


h, (X y) 


The first level of analysis is conducted by setting y = 


~~ we 


y°, a feasible value of y. Then solve 
ae ¥ 50 
minimize: F(X; sy ) 
ea ee ib? 
1 


IV 


subject to: g.(X.,y°) 0 


Ss - 
h; (Xs sy ) = (0, 


9° 


of analysis seeks to find the value of y which produces the 


The solutions are designated Gi and xt The second level 


Mmamimum value of 


ed 


2CG, oe y) = Oe, y) + es OE y) 


Supyect to: CGE > 0 =e 7 
~) ~ 
h,; (Xs ,y) = Or 
Designate the solution by yt. An iterative sequence is 


now established py replacing y~ by ye in the first level 


ww 


problems and resolving; using the resulting Xe to mind. v 
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FIGURE 4. Model Coordination 
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in the second level problem, and so forth until the improve- 
ment achieved with each iteration is less than a specified 
tolerance. 

This approach is called the model coordination method 
because the task of the second level control is to choose the 
linking variables in such a way that the independent first 
level systems are forced to choose solutions which in fact 
correspond to an overall system optimum. In some references, 
this is called the feasibie method. 

The second method, called goal coordination or the dual 
feasible method, views the decomposed system as in Figure 
5. 

It is important to note that in this formulation of the 
problem, . does not necessarily equal a The interactions 
have been literally removed by "cutting" all links between 
subsystems. 

The physical requirement that, in the end, , must equal 
z, termed the interaction-balance principle, 1s satisfied in 
the course of the analysis as follows. 


~e 


In the first level analysis, let i =~°. Then solve 


Pee, . Vek are a5 
minimize: Ly (X, Ty AGA )=£, 0X) +A, T oA, 


subject to: 8, 6X, TT) > 9 
h;(X) Ty.) = 6 
and 
minimize: Ly (X, Ais ‘i n° Eis on =r, ioe + Ne ie 
subject to: £5 (X, Kok. ye 0 
hy (X, XT) = 0 
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ylelding: Ge or and ae 


ow 


The second level problem then becomes choosing A such 
that sclutions to the first level problems result in 
satisfaction of the interaction balance principle. This is a 
well behaved optimization in its own right and is solved with 
the usual techniques of mathematical programming. 

Notice that in this method the coordination effect of 
the second level analysis is effected by manipulating the 
goals of the first level analysis through adjustment of the i 
coordinating variables, hence the term goal coordination 
method. The d multipliers enter the individual first-level 
problem objective functions linearly and act like prices, 
adding to or subtracting from the performance function of each 
subproblem in direct proportion (with proper sign) to the amount 
of Zs demanded and the amount of ve produced. Thus the 
second-level goal coordination can be interpreted as modifying 
"prices" of the interacting variables in order to force the 
independent first-level problems to select consistent values 
of the linking variables and hence the correct overall system 
optimum. 

Much additional information is available in the results 
of the steps of the solution when cast in this format, and the 
interested reader is referred to the considerable literature 
Sueeciemstbject, 113, I¥l915 and 16, for instance]. 

Because of its more straight forward formulation, the 
Sample analysis in the following section is cast in the model 


Soordination format. 
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ie PRELIMINARIES 

The methodoiogy of thermoeconomic analysis can best be 
described through demonstration with a sample analysis. For 
this purpose, an extremely simplified thermal system was 
selected; one which contains the essential features of a 
realistic system but avoids a number of complications which 
would tend to obscure the technique. It should be well under- 
stood that with lumped component representations and several 
Significant losses neglected, the model chosen can not be 
treated as representing a practical plant, nor can the 
results of the analysis be taken as having implications for a 
real system. The model does have many Similarities with ocean 
thermal energy conversion plants as presently conceived, and 
in Appendix B a discussion 1s presented as to what refine- 
ments would be necessary to extend the sample model into one 
of a functional ocean thermal system. For ease of exposition, 
the model will be ciscussed without repeated references to 


these departures from realism. 


B. Daiokemot oa DESCRIPTION 

Figure 6 diagrams the basic system considered. The 
thermal source consists of an infinite supply of seawater at 
a temperature of T = 85°F. The thermal sink is a similarly 


HE 


limitless supply of seawater at Top = uS°F. In the energy 


extraction component the ammonia working fluid is heated as 


Zo 
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FIGURE 6. Basic System 





it flows through the shell side of a rectangular crossflow 
shell and tube heat exchanger with smooth staggered tubes. 
The heat 1s provided by relatively hot seawater flowing 
through the tubes as shown in Figure 7. 

Single phase heat exchange takes place in the device, 
with both fluids remaining compressed liquids. 

The energy conversion component receives hot ammonia 
liquid from the heater, accomplishes energy conversion to 
electrical form, and discharges the liquid at a lower 
temperature. The manner in which the conversion takes 
place is unspecified and not necessary for this sample 
analysis. The conversion process is described by a single 
parameter, Ww, which measures how much energy 1s converted 
to electricity per pound mass of ammonia flowing through the 
device per degree Fahrenheit temperature drop. The value 
selected for W was 0.5 Btu/lbm °F. which is approximately 
half the specific heat for liquid ammonia and, incidentally 
about the same energy available to a turbine with saturated 
vapor inlet conditions.” 

The fluid pressurizer consists simply of one or more 
standard centrifugal pumps, sufficient to drive the working 
fluid through the system at the required rate. Both the hot 


and cold seawater are similarily pumped. 


wy 


“Saturated ammonia vapor at 80°F has enthalpy of 630 
Btu/lbm. Isentropic expansicn to S0°F results in enthalpy 
Peo eee, tam,.or 0.5 Bru/ibm per °F £17]. 
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The system was zoned as depicted in Figure 6, with 
zone 1 consisting of the heat exchanger and working fluid and 
hot seawater pumps, and zone 2 consisting of the energy 
conversion device and the cold seawater pumps. 

Before proceeding further, it should be made clear that 
none of the above assumptions nor those which follow consti- 
tute final arbitrary design selections. Each parameter, 
fluid, and configuration is eventually fixed as an output of 
the analysis itself. Their initial specification should be 
regarded as tentative, pending further information to be 
developed in the course of the study. This preliminary 
configuration acts only as a starting point. 

The next step is to characterize zonal inputs and outputs 
in terms of appropriate physical variables which are descrip- 
tive of the transactions taking place at zone boundaries. 

The principal feature of the hot seawater is its temperature, 
Tip so this was coe as the input to zone 1 from the thermal 
source. The other input to zone 1 is the ammonia discharge 
from zone 2, which is again described by its temperature, i 
The output of zone 1 is hot ammonia liquid at temperature Thy: 
The only remaining variables which cross zone boundaries are 
the electrical output of zone 2, G, and the cold seawater 

from the thermal sink at temperature Tor: 

The global problem is to maximize the profit obtainable 
by selling the system's electrical output at market prices. 
Translated into zone terms, this implies that each zone should 
preduce the required level of output at minimum cost, given 


the inputs it has to work with. 
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ie. Eee) DESCRIPTION OF ZONE 1 

The heat exchanger has length in the direction of sea- 
water flow (2), height in the direction of working fluid flow 
(a), and width transverse to each (w). Tubes have inside 
diameter (d), wall thickness t, and have transverse and 
longitudinal spacing Sp and S;. The inside and outside heat 
flow resistances due to chemical and biological fouling are 
combined into one fouling resistance, Ro. 
Selection of ammonia as the working fluid and seawater 


as the heat source leads to the following table of approx- 


imate physical properties, all considered constant. 


Table IlfI-l. 
PoUtiee ROPER LES 


Sele Joe Ie 
EisuiD DENSITY VES EO CITY SOUVVETEVITY HEAT 
(lbm/ft3) (lbm/ft-hr) (Btu/hr-ft-°F) (Btu/lbm°F) 
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The working fluid mass flowrate (m) and hot seawater mass 
flowrate (my) are provided by centrifugal pumps, which deliver 
the required flows against the head created by frictional and 
form losses in the heat exchanger, (minor losses were 
neglected but could easily be included). The pumping power 


for these pumps, Me and W is a parasitic deduction from the 


Hi 


gross plant electrical output, G.. 
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Fixing the input, output, and linking variables (Tips 


a a TO: and G) temporarily helps to focus an understand- 


ie HH’ 
ime Of zone 1 objectives and constraints. Typ and Top were 
set previously at 85°F and 45°F, and these define the thermal 
potential available to the system. If half of this potential 
is assigned to drive heat through the exchanger surfaces, Ty 
and T. selections of 75°F and 55°F result. If overall power 
output is set at 25 mw (a frequently encountered figure for 


prototype ocean thermal plants), the required working fluid 


flow rate can now be determined from the relationship 
G =m w(T-T). Cal), 


It now becomes evident that the task of zone 1 is to 
receive TuE and qT. and produce the required m at Ty with 
minimum cost. The next major task is to select the design 
variables to be used. To do this, we first look at the 
governing physical and cost relationships. 

As discussed in [18] the performance of a heat exchanger 


mam be deseribed in terms of its effectiveness: 


€ = ne eraser 
(ne C 
where 
ee ae 
and , 
U C 
ye JEE 9 = pH 
mC mC 
p Pp 


Alternatively, the amount of heat transferred can be found from 
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In either case, the fundamental process description is in 
terms of heat exchange surface area, heat exchange coeffic- 
ients, flowrates, temperatures, and fluid properties. 

The major costs in zone 1 are the capital costs of the 


heat exchanger, the pumps, Zn > and the cost of the 


Za 
pumping power. This latter can be considered as an oppor- 
tunity cost and valued at the amount which could have been 


realized had the parasitic pumping power, Me and W been 


H? 
sold at the prevailing rate, Po> instead of being used 
internally. Correlations are available [20] which give 
capital costs of heat exchangers as functions of heat 
exchanger area and capital costs of pumps in terms of the 
product of flowrate and head. The frictional head is usually 
determined empirically and related to flow velocities, 
exchanger configuration, and fluid properties in terms of 
Reynolds numbers. ° 


It might initially appear attractive to choose the 


A and the friction heads 


Uigs Any 


design variables to be Tes 





The Reynolds number is a dimensionless grouping of 
physical variables which indicates the ratio of inertia 
forces to viscous forces. It is formed from the product 
of a characteristic velocity times a characteristic 
dimension, divided by the fluid kinematic viscocity [21]. 


The Prandtl number is also dimensionless, being a 
Measure of the ratio of the diffusivity of momentum to 
Biewdiitusayity of heat. It is formed by multiplying 
the fluid's specific heat times its viscocity and 


ieading by 2ts conductivity (21). 
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DP and DPy> but it is at this point that a key feature of 
thermoeconomic analysis comes into force. Recall that what 
is desired 1s a way to discover the tradeoffs between cost 
and performance. Whatever variables are chosen, it must be 
possible to establish this balance through the functional 
relations over the domain of the variable set. As an 
example of what happens otherwise, consider working with the 
variables suggested just above. Performance can always be 


improved by increasing U and it doesn't cost anything to 


oe 
do so since Ui; is absent from the cost correlation. Costs 
can always be decreased by dropping the DP's, and performance 
would seem to be unaffected because the heat exchange equations 
do not contain pressure drop terms. Any sensible computer 
code would therefore drive U and DP as high and iow, 
respectively, as is allowed. Setting a constraint on these 
parameters is, in effect, arbitrarily choosing them, and no 
information has been gained in the process. Achieving high 
heat transfer with low pressure drops is known to be 
eS eee ec auprion. . 

Besides not permitting a cost-performance balance to be 
weighed, there is a second problem with the variable list 


Suggested, namely that U,, and the DP's are inextricably linked 


H 
through the Reynolds numbers. With a given working fluid, 
heat transfer can only be improved by raising the Reynolds 
numbers with the concurrent result that the pressure drops 


are increased simultaneously. This is the core cencern of 


the zone 1 analysis: UL, is made up of the inside coefficient, 


oy 





h the outside coefficient, ho» the tube wall conductivity 


i? 
and the fouling resistance. There are two pressure drops of 
concern, one in the seawater and one in the working fluid. 
What combination of these parameters will produce the heat 
exchange required at minimum cost? 

To identify an appropriate variable set over which to 
find this cost-performance balance, one must look to the 
next level. Besides physical fluid properties, Reynolds 
numbers depend upon flow rates and spatial dimensions. 

Surface area depends on spatial dimensions. Pump work depends 
on flow rate and pressure drop, which vary as functions of 
Reynolds number. 

Clearly, then, all the cost and performance calculations 
can be built up in terms of flow rates and spatial dimensions, 
and this is the highest levei of variable with which the 
desired tradeoff can be made with the functional relationships 
available. Note, however, that if other valid relationships 
could be Senna which gave the information required in terms of 
other quantities, other variables might be able to be used. 

With the tools at hand, though, it was decided to perform 
the analysis in terms of dimensions and flowrates. But 
variable selection is not yet complete; one has the option of 
which parameters will be allowed to vary independently in the 
code and which will be controlled externally. This question 
is resolved as a matter of judgment, and depends upon the 
confidence the designer has in his preliminary intuition and 


what specific information he seeks from the analysis. Hope- 
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fully, this matter will be clarified as the analysis proceeds. 
For the present investigation, it was decided to permit four 
variables to "float": seawater flow rate, and the length, 
height and width of the heat exchanger. The following 
parameter selections were made. Included are those which 


have been discussed previously. 


Table IIiI-2. 
SEGECTED VALUES OF PARAMETERS 


working fluid: ammonia 

tube diameter: Wie sayin 

tube wall thickness: OR US Se tnch 

tube spacing, Sp? 1.5 d inches 
tube spacing, Sy: 1.5 d inches 
hot seawater temperature: S5°F 

hot working fluid temperature: Hise: 

cold working fluid temperature: Syoaels 

Gross plant power output: 25 MW 

energy conversion factor, pp: Oro bul, bm— = — 
market price of energy: 0.03/Kw-hr 

pump efficiencies: Grd 

fouling heat resistance: 0.005 ft tor Ppt 
tube wall conductivity: 30 Btu/ft-hr-°F 


The functional relationships over both costs and 
performance are developed in detail in Appendix A and summar- 


mZed below. 
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So far, the sample model has been zoned, input and out- 
put variables chosen, design variables selected, and the 
relevant functional relationships formulated. The next step 
is to define explicitly the objective and constraints 
applicable to zone 1 and proceed with computation. The zone 
objective which supports the global objective (maximum profit) 
most directly is acne zone cost. With m fixed by (1), the 
fundamental zone constraint is to raise the required ammonia 


temperature from Le Ou TL given the heat source at Tur By 


H? 


writing this constraint in terms of heat exchanger effective- 


ness, the zone 1 optimization problem may be cast as: 


minimize: BS ae Ly Bi site a Nyy) CA) 
x 
ieee” le 
Stigecet co; & = 2 eee 
HE e 
where 
X = (ms 5 QR, Ws a) 
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The reader may wish to test his own intuition at this 
point. Should the heat exchanger be particularly long, high, 
or wide, or should the dimensions be approximately equal? 

Would seawater flow be expected to be about the same as ammonia 
flow, or more, or less? A very compact heat exchanger would 
reduce the capital costs for that component, give high fluid 
velocities and hence good heat exchange, but would give 

higher pressure drops and hence more pumping costs than a 


larger exchanger. 


D. SOLUTION ALGORITHM 

Numerous computational methods exist for solving (A). 
Since most of the equations are in the Cobb-Douglas form 
f(X) = C ale Caer Pe op ocsuploure use Ceomerrl > 
programming [22] or a number of search techniques. The 
algorithm utilized in this analysis was SUMT4, a complete 
description of which is contained in [23]. A full treatment 
of the underlying theory is given in [12] and only a brief 
explanation of the essential computational sequence is given 
here. 

SUMT4 was chosen partly because of the generality of the 
optimization problem which it can solve. The solution tech- 
niques do not depend on any special features of the problem 
structure, and the conditions for the existence of a solution 
are not overly restrictive. The only absolute requirements 
are that the feasible space be non-empty and that local minima 
Occur at points short of infinity. In order to assure unique- 


ness of the solution, the objective function and all inequality 
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constraints must be continuous and convex/concave respectively, 
but many problems have been solved when these conditions were 
net met [17]. 

This generality is achieved by transforming the con- 
Strained minimization problem into a sequence of unconstrained 
problems, the solutions of which converge to the solution of 
Ehe original problem. The basic idea is that the objective 
and constraints are formed into an auxiliary function 


(penalty function, generalized Lagrangian) as follows: 


given: minimize: f(x) 
Subject to: BX) Ue = le ee. ss Mm 
hs €X) eon) ~ Mis. 5 mt 
where X is an n-dimensional column vector, form: 
Pe ~ m ~ m+p ~ 2 
Beem weet 1 is in BEO.O) o: RS oO. as Ger 
aaa J ies 


Look first at the term involving the inequality 
constraints, B, (x) > 0. Because of the shape of the 
logarithm function, a very large amount is added to the P 
function as the boundary is approached, while the term is 
Ber ateiyeiy flat in the interior of the solution space. 
Conversely, departure in any direction from the equality 
Constraints, Ha ¢ invokes increasing penalties. The size of 
the penalty is controlled by the parameter r. Given a 
weemtine point , Xo Sidwanieinitaal mp, the algorithm searches 
for a minimum of the P function. The search technique is 
specified by the user, and can be either of two modifications 


of the generalized Newton-Raphson method, the method of 
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steepest decent, or McCormick's modification of the Fletcher- 
Powell method. When a minimum has been found to a tolerance 
Specified by the user, the parameter r is reduced by a 
specified ratio and a new minimum is located. This process 
1s repeated until no significant improvement in the P value 
is obtained, as indicated by one of a number of available 
tests. As each successive point is generated, first and 
second order extrapolation are used for convergence 
acceleration. 

Other attractive features of SUMT4 besides its generality, 
have to do with the flexibility it makes available to the 
user. The ability to specify 12 options and 4 tolerances 
afford great freedom in adapting the program to the specific 
application. In addition, the program will produce for 
itself much of the information normally required to be 
supplied. If a feasible starting point is not known, one will 
be calculated. An initial r value may be prescribed, or the 
program will find a good one on its own. If the derivatives 
Sitne £, g, and h functions are not continuous or not 
explicitly available SUMT4 will compute them with central 
differencing procedures. If the problem does have a special 
structure, this can be exploited. A very helpful feature is 
that the mixed interior-exterior penalty function makes it 
unnecessary for the starting point to strictly satisfy the 
equality constraints or be interior to the solution space. 
This makes it possible to avoid much advance manual computa- 
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Finally, SUMT4 is written in a modular structure to 
facilitate changes in logic, options, problems, and input- 
output. This fact was appreciated in the early stages of the 
present investigation, as a few minor modifications weve 
necessary to help avoid programming problems. As an example, 
note that the form of the modeling equations requires that the 
X values be raised to fractional powers. Since Fortran 
accomplishes this through logarithms, negative values of X 
can not be handled. Although SUMT4 will include non-negativity 
constraints automatically, if desired, the differencing and 
extrapolation subroutines still selected negative values during 
the computations. Small changes in the subroutines which 
handle these phases avoided this problem, making the errors 


leading to negative X values easier to locate. 


re INITIAL COMPUTATIONAL CONSIDERATIONS 


The penalty function formed from problem A is: 


~ Ty - oe 
POX.) = eit rly E- =e 
HE eS 


Observe that z represents the costs over the lifetime of the 


plant and is therefore a very large number, while € and 
ye a qT, 
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the algorithm can accomplish major reductions in the P function 


are of the order of one half to one. Without scaling, 


by merely calculating moves which reduce the z values without 
Boe ced ee LOeiew well the equality constraint is satisfied. This 
distortion is easily avoided by scaling the cost equation to 


bring it to approximately the same size as the constraint 
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equation. One preliminary computational run is all that is 
necessary to identify an appropriate scaling factor. 

Secondly, note that the constraint function consists of a 
nested pair of exponentials which is asymptotic to the value 
imo. If the initial x guess produces an € value which is well 
out on the flat portion of the graph, the algorithm is unable 
to determine which direction will produce improvement, since 
the gradient vector is essentially zero. Because of this 
particular shape of the constraint function, it is necessary 
to provide a starting point that is "sufficiently feasible" 
to provide for some gradient in the constraint function 
despite the fact that, in general, SUMT4 does not require 4a 
feasible initial point. This requirement is not particularly 
burdensome, however, because a simple hand calculation can 


owe 


gquaekiy jocate an X, that corresponds with a mid-range e. 


an pinot ZONE 1 SOLUTION 

With the preliminary parameter selections listed in 
iaptoe2 and a starting point of x ovine OemlG. 10). the 
computational sequence was initiated, and the first interesting 
piece of information was developed. With the specified 
weighting of capital costs and costs over lifetime pumping 
charges, the calculation determined that overall costs could 
be continuously reduced by extending the heat exchanger width 
indefinitely. Since this dimension is normal to both fluid 
flow paths, both flow velocities are reduced as w is increased, 
resulting in progressively lower pressure drop losses and 


hence pumping power requirements. The corresponding increase 
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in exchanger capital cost was not sufficient to offset the 
Savings in pump capital cost and power charges. 

This is the first category of information which 
thermoeconomics provides to the designer: which parameters 
should be set as large or as small as other practical consid- 
erations permit. The preferred method of handling instances 
of this sort is to explicitly determine the costs associated 
with allowing the parameter in question to assume increasing 
values. In this case, it should be possible to estimate the 
increase in hull or platform costs required to accommodate 
increasingly wide heat exchangers. If no functional relation- 
ship is available, however, the parameter can simply be 
constrained to some limit which seems reasonable, and the 
Sensitivity of the resulting design to this constraint deter- 
mined. In order to proceed, it was HeCesecadt yar ame Lude a 
wezewlimit in the list of preliminary specifications, so a 
dimensional constraint of 20 feet was chosen, based loosely 
on the arrangement considerations in [16]. 


With this addition, the following solution was obtained: 


My, = 21.1 + 10° 1lbm/hr 
2 = 6.1 ft 

Wea 20.0) ft 

Ble = alee May ae 


This is an example of the second category of information 
obtainable from the method: a point design based on the 
Supporting assumptions and specifications. It is useful to 


look a bit closer at the results and Dipiteadtions Of Enis first 
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run. With the use of a separate small Fortran program, the 
relations in Table III-3 were calculated, using the above 
solution. 

It is important at this stage to insure thet none of 
these variables is out of range with respect to the modeling 
equations. In particular, the Reynolds numbers must be within 
the scope of the friction ioss and heat transfer correlations 
and the elements of the cost equations must lie within the 
span of the cost correlations. If not, the correlations must 
be altered to restore applicability. Comparison with the 
ranges of validity specified in Appendix A shows that all such 
constraints were satisfied in this case. 

Looking first at the physical dimensions which have been 
computed, it is seen that an exchanger has been specified 
which is approximately three tenths as long in the directions 
of flow as in the breadth direction. The immediate implication 
is that the cost tradeoffs specified have dictated that flow 
velocities be reduced, suggesting that pressure losses incur 
greater penalties than the rewards of enhanced heat transfer. 

A possible explanation for this fact can be found through 
an examination of the constituent elements of the overall heat 
transfer coefficient. Note that this number is calculated by 
taking the reciprocal of the sum of the heat flow path 


resistances: 


vy. = aac 2 iL 
H — R. + + RR. + ~ ~GO16+.0013+.005+.00009 
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It is clear that the resistance due to fouling, Rp» is 
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Table I1II-3. 


Puc oweEeoone 2 oOLUL LON 


Heat exchanger length, 2& 
Heat exchanger width, w 
Heat exchanger height. a 


Heat exchanger surface area, Ay 


Number of tubes, Nw 


Heat exchanger capital cost, Zy 
Seawater 


Flow rate, Ma, 

Reynolds number, R OH 
Prandtl number, P oH 
uH 


Heat transfer resistance, R. 


Nusselt number, N 


iieaetion factor, fy 


Friction pressure drop, DP,, 


Friction pump power, WP, 


EMmp Capital cost, Zor 


Working Fluid 


Flow rate, m 
Reynolds number, oo 


PicanGdulemUumMber , a 


Nusselt number, Ny 

Heat transfer resistance, 
BPetettom taetor, f 

Friction pressure drop, DP 
Friction pump power, WP 


Pump capital cost, Zino 


Fouling heat transfer resistance, Re 


Tube wall heat transfer resistance, Ne 


Overall heat transfer coefficient, UL 


ietalezene | cost, z 
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20.0 ft 
6.1 ft 
24,924 ft 
31,368 
$656,711 
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8672 
6.79 
7956 
00162ft--hr-°F/Btu 
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7823 
2.08 
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.00126 ft¢-hr-°F/Btu 
“L2E 
65 1b,/ft* 
5.8 KW 
$341,155 


lbm/hr 
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$1,427,067 





dominating the overall heat transfer coefficient. Even doub- 
ling both of the convective coefficients would only achieve 
an 18 percent improvement in overall transfer. To achieve 
this small improvement, the Reynolds numbers would have to be 
more than doubled, which would increase pumping power 
requirements by almost a factor of 5. 

It should be pointed out that this particular performance 
trade-off constitutes one of the major artificialities of the 
Simplified sample model. In amore realistic system the 
working fluid would be vaporized and the difference between 
the saturation pressures of the hot and cold working fluid 
would represent the bulk of the pump head. Under these 
circumstances, friction losses would not be nearly so control- 
ling, and the heat transfer/pressure loss balance might 
reverse. If the system were as described, however, there is 
an inexorable logic in the dimensional relationships. 

Much additional information can be identified within the 
computed data, some of which will be pointed out later. for 
now, just observe the remarkable balance among the cost 
elements. Not only are the capital costs of the seawater and 
working fluid pumps about the same, but their total is of 
about the same order as the capital cost of the heat exchanger. 
The ratio of three between seawater and working fluid flow 


rates begins to make sense when related to these balanced costs. 


GC. Teo tinGe tHE SOLUTLON 
Wie mete Entiat a solution was Obtained to the initial 


zone problem was most encouraging. 
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The possibility remained, however, that the 
minimum found was not unique and hence did not represent a 
global solution. * Although consideration of the smoothness 
of the functions suggested that the local solution was also 
a global one, further tests were necessary to strengthen that 
conviction. A variety of starting points was accordingly 
devised in order to find out if the same solution resulted. 
Each of the x vectors in Table III-4 led to the identical 
pedwt ion. 

A number of other computationai controls were also 
varied without altering the resulting solution. Tolerances, 
completion criteria, differencing step sizes, scaling factors, 
and minimization methods were specified over wide ranges, and 
the identical solution resulted each time. 

Additionally, SUMT4 contains its own tests for convexity. 
If the matrix of second partial derivatives of the P function 
1s not positive definite at any step in the computation a 
warning message is printed to alert the user that the problem 
is probably not convex (the algorithm proceeds despite this, 
with an orthogenal move). 

This warning message was occasionally received once or 
twice at the beginning of a computational sequence when the 


initial vector was radically different from the eventual 


is this section, the terms local and global are used 
in relation to the zone 1 optimization problem. 
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solution, but any reasonable starting point proceeded 
through the sequence without orthogonal moves being required. 
Finally, belief in the uniqueness of the computed 
solution was supported by the consistent and rational way the 
solution shifted in response to variations in the supporting 
assumptions. In later sections, these variations are 
reported and analyzed. 
Although the above arguments do not constitute a 
Mathematical proof that the prerequisites for solution 
uniqueness exist in the sample model, they do provide strong 


eupport: for such a conclusion. 


lobe m@erntnwt PON PROCEDURES 

Besides the built-in acceleration features of the SUMT4 
algorithm, the user may save substantial calculation time 
with a few simple steps once confidence is gained that the 
local solution is unique. Recall that SUMT4 solves a series 
of minimization problems controlled hy the parameter r. 
Both the starting r and the ratio by which r is reduced in 
each step are under the control of the user. IN Apparel Li ears! 
fairly large initial r (like 1.0) is used when the solution 
1s not known even approximately, and the first trial vector 
may be far from arenes The reduction ratio brings r 
down gradually, with the result that the code can spend 
considerable time solving the wrong problem; i.e., a long 
series of large r subproblems, before arriving in the vicinity 
of the desired solution. Frequently in the course of the 


study however, the analyst wishes to observe the effects of 


oy) 





only small perturbations in the baseline configuration, e.g., 
the response to a 10 percent change in a material price. 
Using the old solution as a starting point, the initial r can 
be chosen small, resulting in only a few minimization sub- 
problems before the new solution is obtained. With a little 
experience, the user can usually select r° and the r reduction 
ratio such that exactly three subproblems are run, yielding 
second-order extrapolation results in the least amount of 
computer time. 

A second control which the user has over computational 
time requirements is exercised through specification of 
subproblem tolerances and completion criteria. These should 
be carefully selected so that computer time is not wasted in 
producing more precision in the solution than is needed. 

With judicious use of the available controls, the 
experienced user can achieve solutions to problems such as 
analyzed here in less than fifteen seconds of computer time. 
Even more oe in computer time may be realized by 
modifying the SUMT4 algorithm for the particular type of 
problem analyzed. Depending on the circumstances, some 
Subroutines may be deleted altogether, others shortened, and 


pamensiton Statements reduced. 


i LINKING VARIABLE BEHAVIOR 

Before proceeding with the zone 1 analysis, and without 
conducting a rigorous second-level analysis as described in 
II C, a measure of preliminary intuition may be gained con- 


cerning the linking variables T,, and T while remaining in the 


H 
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present format. 

Recall that in the first zone 1 solution, the require- 
ment placed on the optimization program was to produce Thy = 
fore dat the minimum cost. This Thy was chosen merely on the 
basis of using half of the available thermal difference for 
heat exchange, evenly divided between the heater and the 
condenser. As far as zone 2 1s concerned, the higher Ty goes 
the better. Raising the enthalpy drop per unit mass of working 
fluid taken across the energy conversion device can do nothing 
but good in zone 2, both from a performance and cost view- 
point. In zone 1, on the other hand, the T,, tradeoffs become 


H 


apparent. To bring T,, closer and closer to the ultimate limit, 


H 
Top? requires increasing heat exchanger effectiveness, either 
through increased exchanger area, better heat transfer 
coefficients, or increased flow rates. Costs attach directly 
or indirectly to each of these improvements, leading to the 
expectation that at zene point it will cost more to raise Ti 
than it is worth in terms of performance. For the first clue 


to optimal T,, selection, therefore, attention was focused in 


H 
zone 1. Leaving all other temporary parameters unchanged, the 
zone 1 analysis was repeated for a range of Ty values. At 
each Ty point, optimum zone 1 X vectors were computed, along 
with zone costs and other intermediate variables. The results 
are given in Table 1, Appendix C. The most important 


relationship developed was that between Ty ele(ol Acigtes IL eereryiee 


shown in Figure 8. 
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The implications of Figure 8 are clear. Certainly, at 
least with the other preliminary parameter assignments, Thy = 


82°F is far superior to Ty = 75°F. As other specifications 
are altered, this will have to be checked, but for the 
present purposes, shifting the Thy specification to 82°F seems 
entirely justified. Without even checking zone 2 fora 
corresponding e, one can also surmise that TO =3 55° bears 


probably too high, anda is of, say, 50°F will be closer to 


the final figure. With these modifications the solution 


became 
M,, = ARTO. lbm/hr 
Miso. Oe tet 
We= 920 it 


and the resulting cost figure was $1,075,991, a reduction of 
25 percent from the baseline case. Rerunning cost against Ty 
at TY = 50°F showed that the minimum cost point occurred at 


Ye 


the same ie as shown in Figure 9. 


we ZONE VARIABLE BEHAVIOR 

All the zone 1 input, output, and linking variables have 
by now been set at least to better figures than arbitrary 
guesses. The next logical step is to refine the interior 
variables. The most arbitrarily chosen of these was d, so 
the next study was concerned with gaining insight over the 


cost effects of tube diameter selection. A series of problems 
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was solved using a range of possible d and t values” with the 
specifications of Table III-2 remaining constant except for 
Ty and TO» which were set at 82°F and 50°F, as discussed 
above. At each d point, the corresponding optimum dimensions 
and seawater flowrates were computed along with the resulting 
costs and other intermediate variables. The results are 
tabulated in Table 2, Appendix C. Again, the most significant 
graph occurs in the cost-diameter plane, as shown in Figure 10. 

In effect, thermoeconomics calls for tube diameters as 
large as possible (another example of type one information), 
requiring the designer to specify a maximum size based on 
practical information not contained in the model. Without 
conjecturing as to what the practical constraint might be, it 
was assumed that good reasons existed for using tubes no 
larger than 2.5 inches. System costs have by now been 
reduced an additional 55 percent, or 66 percent less than the 
baseline case. 

The next parameter studied was Rp- surely lower Re 


values are better than higher ones, so the R, analysis was not 


Ee 
designed to lead to discovery of which way to go. Instead, 
this is the kind of study which develops the third type of 
information recoverable from thermoeconomic analysis: how 
much it is worth to achieve a given technological improvement. 
In other words, 1f an engineer could bring R, from .005 to 


is 


> Tubewall thicknesses were related to tube diameters 
With the guidance of Table 2, Appendix C, reference [19]. 
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.0008 by the addition of special devices or non-fouling 
materials, how much should you be willing to spend for these 
technological advances? The data from the Re runs are given 
in Table 3, Appendix C. Working with these figures leads to 
the following table of total zone 1 costs improvements 
achievable by the reduction of R 


Fr’ 


Table III-5. 
ROW mNGeEACTOR COST EFFECTS 


Be Goon COST SAVINGS* 
.005 U82,566 

004 u72,352 0) al 
.003 461,306 11,046 
.002 Lug 5838 12,468 
.001 u34 592 14,246 
.0009 433,044 1,548 
.0008 u31,470 1,574 
.0007 429,869 1,601 
. 0006 NOMS eye 1,630 
.0005 426,579 1,660 


*Cost Savings for reduction to next lower Rp Resumes 


An Re of 0.0008 was assumed to be all that was techni- 
cally achievable without prohibitive cost and was therefore 
specified in the developing design. This particular value 
also implied that heat transfer would no longer be dominated 


by fouling resistance, since the heat transfer resistance 


values for this latest design were: 


eyll 





Ree OO U0S 


E 
Ry = 0.0004 
Ry =O U0 I 
hep" 0 0023. 
O 


With new d and Re values, it seemed prudent to recheck 


i but, again, the minimum cost occurred at ie =e F anc 

the cost - Th graph had the same shape as in Figures 8 and 9. 
The next variable examined was the arbitrary dimensional 

Bitmiicatton of 20 feet, with the results shown in the following 


table, III-6 and Figure 11. Complete data are given in 


Table 4, Appendix C. 


fabites 1 t=6.. 
DIMENSIONAL CONSTRAINT COST Peeeecrs 
sa eS ee VOL Ur 


TD 


CONSTRAINT COST COST SAVING# 

20 431,470 ase 

25 391,303 40,177 
30 361,406 29,897 
35 337,994 23,412 
40 318,955 18,999 
4S 303,159 15,836 
50 289,686 13,473 


*Cost Saving achievable by extending the constraint by five ft. 


Por the purpose of the Sample analysis, it was assumed that 
the twenty foot constraint remained JL weigh te Te pear 
Fanally, the solution must be tested for SCnS lita Vacay Feo 


the modeling equations themselves. The following table gives 


G2 
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the percent cost change in response to a 10% increase in 


selected parameters. 


Sbevodle IEE Say) = 


COom osENouily TTY TO MODELING EQUATIONS 


PARAMETER 


Py 


Zn 9 equation exponent 


24 equation exponent 


v 
C 
P 


fis equation exponent 


N, equation exponent 


CoH 


i 
NH Sciam xp. 


Z 


f equation exponent 


0 
2O 


f equation constant 


f,, equation constant 


H 


A 
N 


Ny equations const. 


ae 
k 
Po 


k 
W 


equation constant 


equation const. 


Cost 
329,049 
331,850 
335,976 
486,293 
482,772 
390,087 
456,549 
411,139 
WD 
ae el i7 
Lug su70 
ui4,351 
L15,4u9 
417,366 
uu3,359 
u38,715 
438,801 
L2u 278 
Lu US55 
427,934 
434,518 
429,088 
MED DOT 
431,350 
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It 1s interesting to note that three modeling parameters 
produce percentage cost shifts more than double the percentage 
change in the quantities themselves, and that they are all 
fluid properties. Besides suggesting that care be used in 
selecting the proper figures, these results indicate that an 
important refinement would be the use of variable properties 
as fluid state conditions change. 

The least sensitive modeling parameter is tube wall 
conductivity. With the other heat flow path resistances 
being so much higher than that through the wall, this is a 
predictable result. If fouling resistance can be reduced, 
and boiling heat transfer included, however, this may well 
no longer be the case. 

The strong sensitivity to certain fluid properties 
suggested a short side excursion into the cost effects of 
using alternate working fluids. Five of the leading conten- 
ders [37] were analyzed with the results shown in Table III-8. 


Pele eowe 3 
ALTERNATE WORKING FLUIDS 





SEAWATER 

FLOW RATE LENGTH HEIGHT COST 
FLUID 1lbm/hr a ft 8 
Ammonia 3310 a0" 3.1 7.9 431,470 
R-12/31 Mate 6 TO" 7.9 7.3 720,723 
R-500 2068 9 10° 6.8 7.6 680,112 
R-31/114 M8 0 10” 7 ine 692,839 
Propane 605 chet ioe Sat foe Sone of 
Isobutane CTO Che 5.4 jo 591,916 
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Recall, however, that these are strictly zone 1 ccsts, 
and that the working fluid significantly affects the 
performance of the turbine, cycle efficiency, and plant 
arrangements. [24,38] Deen ne stiperiority in zone 1 of 


ammonia as a working fluid is clear. 


Ne ZONE ONE SUMMARY 
At its present stage of development, the system design 


parameters are summarized in Table III-9. 


ihaple wit i—9. 
(PROVE DeEARAMETER SELECTIONS 


Gross plant power output 25 MW 

Working fluid ammonia 

Hot seawater temperature eo Gels 

Cold seawater temperature Ho F 

Hot working fluid temperature als 

Cold working fluid temperature 50 Or 

Tube diameter 2.5 inches 

Tube wall thickness 0.148 inches 

Tube spacing 3.75 inches 

Fouling heat resistance 0.0008 ft*-°F-hr/ 
Btu 

Tube wall conductivity 30 Btu/ft-hr-°F 

Energy conversion factor O- ceeEU/ Lbm=2 EF 

Market price of energy $0.03/kw-hr 


The resulting design is summarized in Table III-10. 


lhe ANALYSIS OF ZONE TWO 
Zone 2 1S analyzed in exactly the same manner as zone l, 


and the information developed is of the same form. Since the 
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Lape  flh—- 10. 


ieNOve econ 1 SOLUTION 


Heat exchanger length, & 

Heat exchanger width, w 

Heat exchanger height, a 

Heat exchanger surface area, Ary 
Number of tubes, Nw 


Heat exchanger capital cost, Z4 


Seawater 


Flow rate, m 


H 
Reynolds number, Roy 
Prandtl number, P 
ielel 


Nusselt number, Nun 

Heat transfer resistance, R. 
Patetlon factor, fy 
Friction pressure drop, DPH 
Friction pump power, Wory 
Pump capital cost, Zany 


Working Fluid 


Flow rate, m 
Reynolds number, ie 
Prandtl number, P 


r 

Nusselt number, Ny 

Heat transfer resistance, Ne 
Beiectton factor, £f 

Friction pressure drop, DP 
Friction pump power, WP 


Rump Capital cost, Za 


Fouling heat transfer resistance, Re 


Tube wall heat transfer resistance, i 


Overall heat transfer coefficient, Us; 


aoecal zone 1 cost, z 
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Gira as ge 


2020 fe 


Peet t 
3298 ft 
1623 

Slew vo 


2 


Soe 10° lbm/hr 


52,648 

aus 

352 
0017 ft°-hr-°F/Btu 
021 

2.1 lbf/ft? 

0.5 kw 

$122,606 


B08 6 We? Gina 


u7,796 
2.08 
278 

MORE 

087 
17.5 lbf/£t? 
1.0 kw 
$116,585 


.0008 ee soe Peren 

0004 ft’-hr-°F/Btu 
192 Btu/ft?-hr-°F 
$431,470 





sample model lacks sufficient realism to make the final actual 
figures meaningful, it was considered redundant to repeat these 
developments pending the intrceduction of more realistic 


two-phase flow conditions into both zones. 


Pl. COMPEETION OF THE ANALYSIS 
With G fixed by the zone 1 analysis, ie determined 


Similarly in zone 2, and preliminary information available as 


memrne range of Re and oe the sample problem could most 
easily be completed through application of the model coordina- 
tion method described in II C. Had the user desired to employ 
the goal coordination method, the zone objective functions 
would have had to be cast in the form required by that approach 
from the outset. In either case, one can see that the 
resulting final design would precisely meet the original 
objective of balancing costs and benefits at the margin 
everywhere in the design. Believing this to be apparent, and 
considering the remaining time available, it was concluded 

that little would be added to the primary goal of this paper 

by presenting the details of the computation. Conclusive 
analysis will be much more valuable if applied to a realistic 
model as described in Appendix B. 

If this is done, a conclusion of great interest will be a 
conclusive determination of an optimum plant output power 
level, G. With increasing returns to scale likely, this 
assignment will probably have to be reached on the basis of 


limitations to the feasibility of commercial manufacture of 
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turbines, platforms, or mooring systems and will accordingly 
not be arrived at trivially. 

Even if this one parameter is selected arbitrarily, 
however, the savings achieved through setting the remaining 
design decision variables at their thermoeconomic optima 


should be considerably in excess of the cost of the analysis. 





TV. SUMMARY AND CONCLUSIONS 


The development in section III does not exhaust the 
possibilities for investigation of the thermoeconomic behavior 
of the sample model. It does, however, indicate the four 


types of information which can be derived: 


Le Which parameters should be fixed as high 
or as low as practical considerations permit. 
Commercial availability or unit costs may 
provide the practical limit, or the constraint 
may be based on technical achievability, as in 
Re» 
arrangements (as in the overall heat exchanger 


or base platform size limitations and hull 


Silvas sLainnage is 


Zz. A point design of the system, based on a given 
specification of parameters which are not 


design variables. 


35 What cost savings could be achieved if the 
practical limits discussed in 1 above are 


extended by a given amount. 
4, The sensitivity of the design to variations 
in the modeling equations themselves. 
Examples of each of these general classes of information 
Ha emeomecimMecdemieseeation Lil K. In addition, Table IV-1 
lists the major thermoeconomics results of the four design 
stages analyzed. 
Certain lessons learned in the process of developing 
these data are considered important enough to warrant emphasis 


through further comment: 
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ise ie Ieee 


SUMMARY OF DESIGN EVALUATION 


VARIABLE a 
m (lbm/hr) 8,500,000 
m,,(1bm/hr) 21,100,000 
merc et) 6.07 
weer t) 2020 
eer to) 3, ins 
Sl am 24,924 
ey (non-dim) BOT? 
eo (non-dim) Wali 28 
U.,(B enre 

migeicu/ ft hr°F) 125 
0, (psf) 65 
Onn (psf) LD 
Z4 Gs) 657,000 
Zn 69) 383,000 
Zno($) 341,000 
Zo ($) 724,000 
Coste (S$) le 27 Olno 


1 - Original design as per table III-2. 


2 - Ty 


Seeede=62.5". ¢ = 0,148" 


= 82°F, T = 50°F 
Le 


4 - R, = 0.098 


F 


2 


5,334,000 

21,964,000 
3.61 
20.0 
6.12 

14,815 

Bao 

9,038 

127 

71 

14 

u74 000 

300,000 

271,000 

571,000 

1,076,000 


aa 


2 


5,334,000 

38,382,000 
3.34 
26.0 
9.83 

4403 

49,086 

bu UBY 

104 

19 

2 

221,000 

130,000 

2e Ono 

253,000 

483,000 


uy 


5,334,000 
33,129,000 
3.11 

20.0 

7.92 
3,297 
52,648 
U7 ,796 
192 
17 
2 
185,000 
123,000 
117,000 
239,000 
431,000 





Realistic, large scale systems (such as discussed in 
Appendix B) are undoubtedly best handled through a zone 
approach. Analytical tools are available to insure that the 
solutions to the zone sub-problems are coordinated such that 
an overall system optimum is achieved. Vapor power systems 
lend themselves readily to the zone approach in that the 
functions of the major components are distinct and the linking 
relationships are clear. 

There is a simple rule for deciding which zone variables 
should be free to vary in the optimizing algorithm and which 
should be externally controlled. If the cost effects of 
varying a parameter can be readily included in the statement 
of the objective function, the parameter may be included in 
the vector of decision variables. If not, the parameter 
should be fixed during the minimization search and investi- 
gated separately. Sensitivity analysis reveals which of these 
affect the design sufficiently that they require careful 
selection. Those that do should be examined further, either 
to develop approximate costing relationships or to find other 
valid means for establishing their final specification. 

The SUMT4 algorithm was found to be convenient and 
effective in finding solutions to the zone optimization 
problems. Although not demonstrated, it is expected that it 
would be equally capable of handling the second level 
coordination problem. Used as a subroutine, it should be 
straightforward to mechanize SUMT into a master program for 


conducting the iterative first level-second level procedure 
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discussed in section III E. Substantial economies in 
computation time appear possible of achievement by adapting 
the general program to the specific application. 

The principal conclusion of this study is that 
thermoeconomic analysis of the kind suggested is capable of 
producing insights which are of considerable value to the 
system designer. This is particularly true when the system 
in question is such that the long run cost effects of technical 
design specifications are not well known on the basis of 
extensive previous experience with similar systems. Environ- 
mental power systems are such a class, and the financial 
imperatives which result from their considerable size, coupled 
with the recognized need for their early success, strongly 
Suggest the wisdom of carefully applying the integrated 


cost/performance analysis techniques presented in this paper. 
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APPENDIX A 
Dey bOrMENr Of fie THERMOECONOMIC MODEL 


1tZ FLUID FLOW AND HEAT EXCHANGE EQUATIONS 


A. Basic Model Description 

The basic model is described in sections II B and C 
of the text. Ammonia was chosen as the working fluid because 
of its favorable heat exchange properties [24], but also 
note the practical objections to ammonia contained in [25]. 
From among the candidate thermal sources and sinks, those 
involved in the ocean thermal energy conversion process [26] 
provide the low thermal potential typical of environmental 


energy sources and were therefore chosen for the sampie model. 


B. Pump Work, Working Fluid (See Knudsen and Katz [27]) 
The work required to pump the working fluid against 
friction losses in the heat exchanger is 


Wo aad ft - 1b, 
p e) ni; el 


where DP depends on a friction factor which is measured 
empirically for the given configuration: 
2 out 
ie Cy N 1b, 


Deo a 
p Bo pe 


For the specific application considered, 


oe. ot ibm ___ 

H W _-d 3600 eee 

Nee = = dimensionless 
L 
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and 


2 Se res Pee dimensionless 
G0 G < Rees ZS O00 ) 
where , 
m Sep 
R 2 =e dimensionless 
ec ULw 
c.. Pump Work, Hot Seawater (See Streeter [28]) 


The pump work required to pump the hot seawater 


against friction losses thrceugh the heat exchanger tubes is 


_— My DP, 
Ho Py Ne 
we th 9 
Gree. 1G 
De =f _ Hs. 
H H 2d 
© ik, Boe 
U_ = HTL 
fF 0,4 wa 
and 
-1/4 5 
oe = R R 
fy Ol, SOLE oH ( oH < 10°) 
where ; 
R Pee Tk 
eH T Up dwa : 
Be Heat Exchange Equations 
The overall heat exchange coefficient is given by: 
_ ae See 
ee Peter Ff Ro + R 
a O W F 
where: 
ilies R.. is the thermal resistance due to fouling 


E 


and is an externally controlled parameter in the analysis. 
Di RR is the thermal resistance due to the 


conduction path through the tube walls, and is computed from 


is 





ees 
W ai 
mw w 
oe R. is the thermal resistance at the seawater/tube 


interface, =. The immer film coefficient, h; is given by 





ai 
[29] 
Inl are | 
tf - z 0.8 ies d<0.055 
ake - Nu, yield Sus R oH Pri (7) 
where . 
um, SAS 
4 Isl Pag 
Roy — TW, dwa_ a US) before, 
and 
C u 
Pr. = —PH HH 
tal k 
H 
TR Ro 1s the thermal resistance at the ammonia/tube 
interface, Ay . The outer film coefficient, hos 1s given 
A 
by [29] ee 
ih d : 
26 = 0.512 R 0.5625, lye 
k ec 
where i 
m Sy 
Roc 2 Tee before, 
and 
C_u 
Pr - _P_ 
k 


mis COST EQUATIONS 

The costs accounted for in sample model were the capital 
cost of the heat exchanger and pumps and the opportunity cost 
of the parasitic pumping power. Rather than encumber the 
model with a detailed accounting of discounted cash flows, ea 


Simple financing method was assumed which was believed to 


76 





provide cost assignments which were adequate for demonstration 
purposes. Specifically, amortization and interest expense 
were assigned as if the capital investment were funded with 
a twenty year, ten percent mortgage. The twenty year life 
was based on the consensus of the working group for economics, 
ocean thermal energy conversion workshop, held under the 
auspices of the National Science Foundation in Washington, 
D.C. in September, 1974 [30]. The ten percent interest rate 
corresponds to the discount rate specified for all public 
programs by current government directives. All costs were 
expressed in terms of the total cost over the twenty year 
life of the plant. Capital costs were estimated from the 
data contained in [20] and multiplied by a factor of 2.3 to 
approximate current dollars. 
A. Capital Cost, Heat Exchanger 

The base cost (BC) graph for shell and tube heat 
exchangers given in the reference was converted into 
functional form as 


OF 627 
2 _ 1d Law 
Be. 1a) Ais 5 where As = 8,5; - 


Applying appropriate adjustment factors in accordance with 
the instructions contained in the reference, the cost 
estimating relationship chosen to represent the capital cost 


of the heat exchanger became 


Weo2/ 


Exchanger Cost = Z, = 500 As : 


Z 





BE. Capital Cost, Pumps 
The base cost graph for centrifugal pumps and 
drivers translates into 


OA (s0) 
BC = 84 (C*H) 


Applying material adjustment factors for monel (seawater 


pumps) and stainless steel (ammonia pumps) the costs become 


Zo = Zo, * Zoo 

; 0.602 
Zo, = 271 (CH,) 

; 0.602 
zee) (CH) 


where the CH terms are the product of the flow rate in gallons 


per minute and the pump head in pounds per square inch. 


C. Cost of Pumping Power 


Since the total system cost was formulated as 


COST = p(W + Wy) + 2) + 2, 


it was necessary to use a Py factor which brought the power 
cost term into the same units as the z terms; 1.e., cost over 
a twenty year period. Applying the necessary unit conversion 


factors, P, was used as 


$ - hr 
= Q Ae eee 
Bg 0.0019795 —— Tb, 


ZUR yrs. 
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APPENDIX B 
EXTENSION TO A REALISTIC SYSTEM 


The sample analysis of section III indicates that valuable 
information could be obtained if the thermoeconomic method 
of anlysis were applied to a realistic system model. This 
Appendix outlines suggestions for modifying and extending the 
sample model so that the insights developed will be meaning- 
ful in the design of an actual working prototype. 

Several groups have already conducted extensive study of 
the ocean thermal energy conversion process, leading to a 
number of preliminary designs [31, 32]. Probably the most 
complete design from the viewpoint of detailed system 
modeling is that of the University of Massachusetts (Amherst) 
[31]. With the engineering relationships already developed 
mene pPOIint of a comprehensive, coordinated, feasible design, 
this model is a very attractive candidate for thermoeconomic 


analysis. 


A. ENGINEERING MODIFICATIONS 

The most dramatic modification to the sample model results 
when two-phase heat transfer is included. The difference in 
Seuubation pressures at ie ess 2h215 TehoYal T. = 50°F amounts to 
almost 10,000 lbf/ft?, dwarfing the 17 Lbf/ft? EOTet ion 
pressure drop in the sample model. The obvious result is that 


reduction in flow velocity will no longer be the dominant 


objective of the optimization search and it may well result 
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that higher velocities are more advantageous than lower ones. 

The second major change resulting from the inclusion of 
two-phase flow is in the area of heat exchange calculations. 
It 1s well known that boiling heat transfer coefficients are 
much higher than those for compressed liquids. With the 
working fluid heat transfer resistance becoming negligible 
in comparison with the other path resistances, further major 
alterations in the optimum design can be anticipated. Tube 
wall resistance may become a significant contributer and will 
have to be carefully specified. 

Again, much of the background engineering in this area 
has already been done. The Amherst team has adapted the Chen 
[33] and Chawla [£34] correlations for boiling heat transfer 
to the plate-fin exchanger configuration they recommend. 

A second refinement which may strongly influence the 
thermoeconomic balance is the inclusion of the multiple 
fluid path flow losses commonly referred to as minor losses. 
Preliminary calculations show that these may well substantially 
exceed the friction and form losses which were the only flow 
losses accounted for in the sample model. Added to the 
substantial power required to pump the cold seawater from a 
depth of some 2000-3000 ft, pumping costs are certain to 
shift the zone costing relationships radically. 

The sensitivity of the sample model to the physical 
characteristics of the working fluid suggests that greater 
attention be given to its selection. As reported by [35], 
these properties also have a strong influence on turbine 


design. 
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B'. COSitING MODIFICATIONS 

First, the tube and shell exchanger correlation used in 
the sample analysis will have to be modified for applicability 
to the plate-fin configuration. Second, note that the pump 
cost estimating relationship employed was for centrifugal 
pumps in ordinary industrial use. The pumps in an ocean 
thermal plant, which are required to drive very high volume 
rates of flow against very small pump heads, will necessarily 
be quite different. Something on the order of shrouded 
propellors in axial flow will be better suited to the 
application, and these will have different cost estimating 
formulas altogether. 

In both cases, much greater care will be required in 
selecting appropriate materials and insuring commonly valued 
dollars than was taken in the model of section III. 

But probably the most difficult costing relationships 
to estimate will be those associated with the size of the hull 
or platform required to support the plant. These costs control 
the component dimensions which can be chosen, and will probably 
dictate the very important ultimate selection of optimum 


plant power output. 


c. FINANCING CONSIDERATIONS 

Rather than assuming a single financing method, the 
realistic analysis should apply the general methods for 
evaluating investment decisions. The overall cash flow 
inherent in the project should be estimated as to its time 


dimension and discounted to its net present value. The 


8}. 





assumptions recommended by [36] should be utilized in per- 


forming this step. 


Be COMPUTATIONAL CONSIDERATIONS 

It 1S recommended that any realistic analysis be de- 
composed into at least three zones grouped about the three 
major functional components: the boiler, turbine, and 
condenser. The linking variables should consist of the 
working fluid state variables at the entrance and exit to 
these components. In addition to the temperatures, it will 
probably be necessary to include vapor quality as a linking 
variable, in that the performance of all three components 
is affected by the moisture content. 

Either the model or goal method may be used in 
establishing problem coordination. Model coordination is th 
more straightforward of the two, but the auxiliary information 
produced in the process of goal coordination should be 
considered in making the selection. 

In either case, the iterative procedure required should 
be mechanized into a controlling program, with SUMT or some 
other appropriate optimization algorithm being applied 


sequentially to the first and second level problems. 
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